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Abstract 


This  Technical  Memorandum  presents  an  evaluation  baseline  for  ship  detection  employing 
RADARS AT-2  Fine  Quad  (FQ)  Mode  imagery  and  an  initial  investigation  of  the  potential 
performance  for  data  obtained  by  a  compact  polarimetry  (CP)  SAR  system.  A  CP  SAR  system 
with  circular  polarization  on  transmission  and  two  orthogonal  linear  polarizations  on  receive  was 
simulated  using  RADARSAT-2  FQ  data.  Polarimetric  SAR  (PolSAR)  ship  detection  algorithms 
were  applied  to  both  the  FQ  and  simulated  CP  data.  From  statistical  decision  theory,  the 
likelihood  ratio  test  with  Neyman-Pearson  criterion  was  used  to  define  a  decision  variable. 

FQ,  dual  polarization  (including  simulated  CP,  and  horizontal  polarization  on  transmit  and 
horizontal  and  vertical  polarizations  on  receive  (FTH  +  F1V)  with  amplitude  only),  and  single 
polarization  (HH)  images  of  known  ships  were  considered.  In  this  study,  the  smallest  detected 
ship  had  a  length  of  1 8  m.  The  results  demonstrate  that  a  quad  polarization  system  provides  the 
best  ship  detection  performance  in  this  application  as  compared  to  single  and  dual  polarization 
SAR  systems.  The  results  also  demonstrate  that  a  CP  system  provides  better  performance  than  a 
conventional  dual  polarization  system  and  in  turn  than  a  single  polarization  system.  It  is 
important  to  note  that  CP,  dual  polarization  and  single  polarization  systems  have  a  swath  width 
that  is  twice  as  wide  as  that  of  a  PolSAR  having  the  same  resolution. 

Detection  performance  was  characterized  by  calculating  receiver  operating  characteristics  for 
RADARSAT-2  FQ  data  and  simulated  CP  data  to  quantify  the  trade-off  between  the  probability 
of  missed  detection  and  the  probability  of  false  alarm.  SAR  detections  were  validated  by  using 
automatic  identification  system  (AIS)  data. 

In  further  studies,  detection  algorithms  based  on  the  likelihood  ratio  test  have  been  optimized  and 
the  decision  variable  defined  using  the  likelihood  ratio  test  has  been  input  into  a  K-distribution 
constant  false  alarm  rate  (KCFAR)  detector  as  the  required  image.  The  KCFAR  PolSAR 
(KCFAR  Quad)  detector  found  all  targets  for  which  AIS  data  were  available,  and  also  additional 
targets.  Results  show  that  this  detector  is  capable  of  detecting  smaller  targets  than  those  for  a 
single  polarization  system. 
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Resume 


Le  present  memoire  technique  porte  sur  une  reference  d’ evaluation  pour  la  detection  des  navires 
au  moyen  d’imagerie  en  mode  a  polarisation  en  quadrature  (quad-pol)  fin  de  RADARSAT-2  ainsi 
que  sur  une  premiere  etude  du  rendement  potentiel  d’un  systeme  de  radar  a  synthese  d’ouverture 
(SAR)  a  polarimetrie  compacte  (CP).  Un  systeme  de  SAR  CP  ou  remission  se  fait  avec  une 
polarisation  circulaire  et  la  reception  avec  deux  polarisations  rectilignes  orthogonales  a  ete  simule 
au  moyen  de  donnees  du  mode  polarisation  en  quadrature  fin  (FQ)  de  RADARSAT-2.  Des 
algorithmes  de  detection  de  navires  par  SAR  polarimetrique  (PolSAR)  ont  ete  appliques  aux 
donnees  FQ  et  aux  donnees  CP  simulees.  Le  test  de  rapport  de  vraisemblance  avec  critere  de 
Neyman-Pearson,  issu  du  de  la  theorie  statistique  de  la  decision,  a  servi  a  definir  une  variable  de 
decision. 

L’etude  a  porte  sur  des  images  de  navires  connus  obtenues  en  mode  FQ,  en  mode  polarisation 
double  (y  compris  en  mode  CP  simulee  et  en  mode  a  emission  en  polarisation  horizontale  et  a 
reception  en  polarisation  verticale  (HH  +  HV)  avec  amplitude  seulement)  et  en  mode  polarisation 
simple  (HH).  Dans  la  presente  etude,  le  plus  petit  navire  detecte  avait  une  longueur  de  18  m.  Les 
resultats  montrent  qu’un  systeme  a  polarisation  en  quadrature  offfe  le  meilleur  rendement  de 
detection  des  navires  dans  cette  application  que  les  systemes  SAR  a  polarisation  simple  ou 
double.  Les  resultats  demontrent  egalement  qu’un  systeme  CP  permet  d’obtenir  un  meilleur 
rendement  qu’un  systeme  a  polarisation  double  classique,  qui  offfe  a  son  tour  un  rendement 
superieur  a  celui  d’un  systeme  a  polarisation  simple.  II  est  important  de  noter  que  les  systemes 
CP  et  les  systemes  a  polarisation  simple  ou  double  ont  une  bande  au  sol  deux  fois  plus  large 
qu’un  systeme  PolSAR  de  meme  resolution. 

Le  rendement  de  detection  a  ete  caracterise  en  calculant  la  courbe  ROC  ( Receiver  Operating 
Characteristics )  a  partir  des  donnees  FQ  de  RADARSAT-2  et  des  donnees  CP  simulees  afin  de 
quantifier  le  compromis  entre  la  probabilite  de  detection  manquee  et  la  probabilite  de  fausse 
alarme.  Les  detections  SAR  ont  ete  validees  au  moyen  des  donnees  du  systeme  d'identification 
automatique  (A1S). 

Au  cours  d’ etudes  subsequentes,  des  algorithmes  de  detection  fondes  sur  le  test  de  rapport  de 
vraisemblances  ont  ete  optimises,  et  la  variable  de  decision  definie  au  moyen  du  test  de  rapport  de 
vraisemblances  a  ete  entree  dans  un  detecteur  a  taux  de  fausse  alarme  constant  a  distribution  K 
(KCFAR)  en  tant  qu’image  requise.  Le  detecteur  KCFAR  PolSAR  (KCFAR  Quad)  a  detecte 
toutes  les  cibles  pour  lesquelles  des  donnees  AIS  etaient  disponibles,  ainsi  que  d’autres  cibles. 
Les  resultats  montrent  que  ce  detecteur  peut  detecter  des  cibles  beaucoup  plus  petites  qu’un 
systeme  a  polarisation  simple. 
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Executive  summary 


Ship  detection  using  RADARSAT-2  Fine  Quad  Mode  and 
simulated  compact  polarimetry  data 

Chen  Liu;  Paris  W.  Vachon;  Ryan  A.  English;  Nicholas  Sandirasegaram 
DRDC  Ottawa  TM  2009-285;  Defence  R&D  Canada  -  Ottawa;  February  2010. 

Introduction:  RADARSAT-2  capabilities  include  polarimetric  synthetic  aperture  radar  (PolSAR) 
modes.  Improved  ship  detection  using  PolSAR  data  over  single  channel  data  has  been  previously 
demonstrated  using  airborne  PolSAR  data.  In  this  Technical  Memorandum,  the  application  of 
PolSAR  ship  detection  algorithms  developed  using  airborne  data  to  RADARSAT-2  Fine  Quad 
(FQ)  Mode  data  is  evaluated. 

Compact  polarimetric  (CP)  SAR  is  an  approach  to  dual  polarization  SAR.  One  form  is  hybrid- 
polarity  in  which  a  circularly  polarized  wave  is  transmitted  and  two  orthogonal,  mutually- 
coherent  linearly  polarized  waves  are  received.  CP  SAR  has  a  wider  swath  width  than  quad 
polarization  systems  of  the  same  resolution.  As  described  in  this  report,  hybrid-polarity  SAR 
using  RADARSAT-2  FQ  data  is  simulated.  PolSAR  ship  detection  algorithms  then  were  applied 
to  the  simulated  CP  SAR  data. 

The  dual  polarization  HH  +  HV  with  amplitude  only  and  single  polarization  (HH)  data  were 
extracted  from  the  RADARSAT-2  FQ  data. 

Based  on  statistical  decision  theory,  the  likelihood  ratio  test  with  Neyman-Pearson  criterion  was 
used  to  define  a  decision  variable.  Detection  performance  was  characterized  by  calculating 
receiver  operating  characteristics  for  RADARSAT-2  FQ  data  and  simulated  CP  data  to  quantify 
the  trade-off  between  the  probability  of  missed  detection  and  the  probability  of  false  alarm. 

The  RADARSAT-2  FQ  data  acquired  over  the  Strait  of  Gibraltar  are  accompanied  by 
simultaneously  acquired  automatic  identification  system  (AIS)  data.  The  targets  detected  in  quad 
polarization  and  simulated  CP  SAR  data  were  validated  using  the  AIS  data. 

Quad  polarization,  dual  polarization  (including  simulated  CP  and  HH  +  HV  with  amplitude  only), 
and  single  polarization  (HH)  images  of  known  ships  were  then  considered:  the  relative  detection 
performance  of  these  various  SAR  systems  was  estimated  by  determining  the  receiver  operating 
characteristics  (ROC). 

In  further  work,  detection  algorithms  based  on  a  likelihood  ratio  test  have  been  optimized  and  the 
decision  variable  derived  using  the  likelihood  ratio  test  has  been  sent  to  a  K-distribution  constant 
false  alarm  rate  (KCFAR)  detector  to  perform  the  pixel-based  detection.  The  algorithms  for 
calculating  the  decision  variable  from  quad  polarization  data  was  implemented  by  the  Defence 
R&D  Canada  -  Ottawa  (DRDC  Ottawa)  using  the  Analysts’  Detection  Support  System  (ADSS), 
which  was  developed  by  the  Defence  Science  and  Technology  Organization  (DSTO,  Australia). 
The  KCFAR  detector  used  in  this  study  was  developed  by  the  Defence  Science  and  Technology 
Laboratory  (Dstl,  U.K.)  in  ADSS.  Initial  results  demonstrate  that  KCFAR  PolSAR 
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(KCFARQuad)  detection  algorithms  can  potentially  be  effective  for  detection  of  small  ships  (1  ~ 
18  m). 

The  KCFAR  Quad  detection  algorithms  described  in  this  report  provide  a  performance  baseline 
for  PolSAR  ship  detection. 

Results:  The  results  clearly  demonstrate  that  significant  improvement  in  ship  detection  is 
possible  when  using  RADARSAT-2  FQ  data  as  compared  to  dual  polarization  and  single 
polarization  modes.  The  simulated  CP  SAR  data  provide  significant  detection  improvement  over 
single  polarization  systems,  as  well  as  conventional  dual  polarization  with  amplitude  only.  The 
detection  performances  of  quad  polarization,  dual  polarization  and  single  HH  polarization  have 
been  estimated.  A  total  of  eight  ships  in  one  scene  were  studied,  for  a  probability  of  false  alarm 
of  10‘4;  the  quad-pol  system  achieved  seven  times  lower  probability  of  missed  detection  than 
single  polarization  systems;  the  CP  system  achieved  four  times  lower  probability  of  missed 
detection  than  single  polarization  systems.  For  the  CP  system,  the  detection  performance  for  five 
of  the  ships  was  significantly  better  than  that  of  the  dual  polarization  system  with  amplitude  only. 
For  the  other  three  ships,  the  performance  of  the  dual  polarization  system  was  slightly  better. 
Overall,  the  detection  performance  of  the  CP  SAR  system  was  judged  to  be  better  than  that  of  the 
dual  polarization  system  while  keeping  the  same  swath  width. 

The  results  demonstrate  the  significant  performance  benefit  of  the  CP  SAR  as  compared  to  the 
single  polarization  HH  system,  while  providing  a  wider  swath  coverage  than  the  quad 
polarization  system. 

For  the  cases  studied,  the  smallest  detected  ship  had  a  length  of  18  m  and  a  speed  of  0.59  m/s. 

The  KCFAR  Quad  detector  provides  good  results.  Thirty  ships  verified  with  A1S  data  in  the 
ocean  region  in  the  Strait  of  Gibraltar  image  were  detected. 

Significance:  The  PolSAR  ship  detection  algorithms  using  amplitude  and  phase  information  of  a 
target  in  four  channels  show  a  significant  improvement  in  ship  detection.  A  unique  feature  of  the 
phase  information  of  targets  obtained  from  a  PolSAR  system  has  been  used  to  reduce  the  false 
alarm  rate  by  discriminating  between  images  of  the  ships  and  azimuth  ambiguities.  The  PolSAR 
system  can  be  used  to  detect  smaller  ships  than  dual  polarization  or  single  polarization  systems. 

The  RADARSAT  Constellation  Mission  (RCM)  will  provide  CP  as  an  operational  mode,  which 
could  be  beneficial  to  ship  detection  activities.  It  is  recommended  that  the  CP  mode  be  considered 
for  wide  area  surveillance,  in  particular,  for  ship  detection. 
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Ship  detection  using  RADARSAT-2  Fine  Quad  Mode  and 
simulated  compact  polarimetry  data 

Chen  Liu;  Paris  W.  Vachon;  Ryan  A.  English;  Nicholas  Sandirasegaram 

DRDC  Ottawa  TM  2009-285;  R  &  D  pour  la  defense  Canada  -  Ottawa;  Fevrier 
2010. 

Introduction  :  RADARSAT-2  permet  notamment  d’utiliser  des  modes  de  radar  a  synthese 
d’ouverture  polarimetrique  (PolSAR).  Des  donnees  de  PolSAR  aeroporte  ont  permis  de 
demontrer  que  ce  systeme  offre  un  meilleur  rendement  de  detection  des  navires  qu’un  systeme  a 
un  seul  canal.  Dans  le  present  memoire  technique,  on  a  evalue  Tapplication  a  des  donnees  du 
mode  a  polarisation  en  quadrature  fin  (FQ)  de  RADARSAT-2  d’algorithmes  de  detection  de 
navires  par  PolSAR  qui  ont  ete  mis  au  point  au  moyen  de  donnees  aeroportees. 

Le  SAR  a  polarimetrie  compacte  (CP)  est  une  approche  de  SAR  a  polarisation  double.  Une 
methode  possible  de  CP  est  l’utilisation  d’une  polarisation  hybride,  ce  qui  consiste  a  emettre  des 
ondes  en  polarisation  circulaire  et  a  recevoir  deux  ondes  mutuellement  coherentes  avec  deux 
polarisations  rectilignes  orthogonales.  Un  SAR  CP  a  une  bande  au  sol  plus  large  qu’un  systeme  a 
polarisation  en  quadrature  de  meme  resolution.  Dans  le  present  rapport,  le  SAR  a  polarisation 
hybride  est  simule  au  moyen  des  donnees  FQ  de  RADARSAT-2.  Des  algorithmes  de  detection 
des  navires  par  PolSAR  sont  ensuite  appliques  aux  donnees  de  SAR  CP  simulees. 

Les  donnees  en  polarisation  double  (HH  +  LTV)  avec  amplitude  seulement  et  en  polarisation 
simple  ont  ete  extraites  des  donnees  FQ  de  RADARSAT-2. 

Le  test  de  rapport  de  vraisemblance  avec  critere  de  Neyman-Pearson,  issu  de  la  theorie  statistique 
de  la  decision,  a  servi  a  definir  une  variable  de  decision.  Le  rendement  de  detection  a  ete 
caracterise  en  calculant  la  courbe  ROC  ( Receiver  Operating  Characteristics )  a  partir  des  donnees 
FQ  de  RADARSAT-2  et  des  donnees  CP  simulees  afm  de  quantifier  le  compromis  entre  la 
probabilite  de  detection  manquee  et  la  probabilite  de  fausse  alarme. 

Les  donnees  FQ  de  RADARSAT-2  acquises  au-dessus  du  detroit  de  Gibraltar  sont  accompagnees 
de  donnees  du  systeme  d'identification  automatique  (AIS)  acquises  en  meme  temps.  Les 
detections  de  cibles  realisees  au  moyen  des  donnees  a  polarisation  en  quadrature  et  des  donnees 
SAR  CP  simulees  ont  ete  validees  a  Taide  des  donnees  AIS. 

L’etude  a  porte  sur  des  images  de  navires  connus  obtenues  en  mode  FQ,  en  mode  polarisation 
double  (y  compris  en  mode  CP  simulee  et  en  mode  a  emission  en  polarisation  horizontale  et  a 
reception  en  polarisation  verticale  (HH  +  HV)  avec  amplitude  seulement)  et  en  mode  polarisation 
simple  (HH).  Le  rendement  de  detection  relatif  des  differents  systemes  de  SAR  a  ete  evalue  en 
determinant  la  courbe  ROC  ( Receiver  Operating  Characteristics). 

Au  cours  d’etudes  subsequentes,  des  algorithmes  de  detection  fondes  sur  le  test  de  rapport  de 
vraisemblances  ont  ete  optimises,  et  la  variable  de  decision  definie  au  moyen  du  test  de  rapport  de 
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vraisemblances  a  ete  entree  dans  un  detecteur  a  taux  de  fausses  alarmes  constant  a  distribution  K 
(KCFAR)  en  vue  d’effectuer  la  detection  simple  par  pixels.  Les  algorithmes  de  calcul  de  la 
variable  de  decision  a  partir  des  donnees  de  quadruple  polarisation  ont  ete  appliques  par  R  &  D 
pour  la  defense  Canada  -  Ottawa  (RDDC  Ottawa)  au  moyen  de  l’Analysts’  Detection  Support 
System  (ADSS),  mis  au  point  par  la  Defence  Science  and  Technology  Organization  (DSTO, 
Australie).  Le  detecteur  KCFAR  utilise  dans  le  cadre  de  la  presente  etude  a  ete  developpe  par  le 
Defence  Science  and  Technology  Laboratory  (DSTL,  R.-U.)  dans  TADSS.  Les  premiers  resultats 
montrent  que  les  algorithmes  de  detection  PolSAR  KCFAR  (KCFARQuad)  pourraient  etre 
efficaces  pour  la  detection  de  petits  navires  (1®18  m). 

Les  algorithmes  de  PolSAR  KCFAR  decrits  dans  le  present  rapport  serviront  de  reference  de 
rendement  pour  la  detection  de  navires  par  PolSAR. 

Resultats  :  Les  resultats  montrent  clairement  que  l’utilisation  de  donnees  FQ  de  RADARSAT-2 
permet  d’ameliorer  grandement  la  detection  des  navires  par  rapport  aux  modes  a  polarisation 
double  ou  simple.  Les  donnees  de  SAR  CP  offrent  une  amelioration  marquee  du  rendement  de 
detection  par  rapport  aux  systemes  a  polarisation  simple  et  aux  systemes  a  polarisation  double 
classiques  avec  amplitude  seulement.  Le  rendement  de  detection  des  modes  a  polarisation  en 
quadrature,  a  polarisation  double  et  a  polarisation  simple  LLH  a  ete  estime.  La  scene  etudiee 
comportait  au  total  8  navires.  Pour  un  taux  de  fausse  alarme  de  1 0“4,  le  systeme  a  polarisation  en 
quadrature  a  obtenu  une  probabilite  de  detection  manquee  7  fois  plus  basse  que  le  systeme  a 
polarisation  simple,  et  le  systeme  CP  a  obtenu  une  probabilite  de  detection  manquee  4  fois  plus 
basse  que  le  systeme  a  polarisation  simple.  Le  systeme  CP  a  obtenu  un  rendement  bien  superieur 
a  celui  du  systeme  a  polarisation  double  avec  amplitude  seulement  pour  5  des  8  navires.  Pour  les 
trois  autres  navires,  le  rendement  du  systeme  a  polarisation  double  etait  legerement  superieur. 
Dans  l’ensemble,  le  systeme  de  SAR  CP  offre  un  rendement  juge  meilleur,  tout  en  gardant  la 
meme  largeur  de  bande  au  sol. 

Les  resultats  demontrent  que  le  SAR  CP  possede  un  avantage  de  rendement  marque  sur  le 
systeme  a  polarisation  simple  1111,  tout  en  foumissant  une  bande  au  sol  plus  large  que  le  systeme 
a  polarisation  quadratique. 

Pour  les  cas  etudies,  le  plus  petit  navire  detecte  avait  une  longueur  de  18  m  et  une  vitesse  de 
0,59  m/s. 

L’algorithme  KCFAR  Quad  donne  de  bons  resultats  de  detection  :  dans  l’image  de  la  region  du 
detroit  de  Gibraltar,  33  navires  dont  la  presence  dans  le  secteur  a  ete  verifiee  au  moyen  de 
donnees  A1S  ont  ete  detectes. 

Portee  :  Les  algorithmes  de  detection  de  navires  par  PolSAR  utilisent  l’information  d’amplitude 
et  de  phase  d’une  cible  provenant  de  quatre  canaux  PolSAR  afm  d’obtenir  une  amelioration 
marquee  de  la  detection  des  navires.  Une  caracteristique  unique  de  l’information  de  phase  des 
cibles  obtenue  d’un  systeme  PolSAR  a  servi  a  reduire  le  taux  de  fausse  alarme  en  eliminant  les 
ambigui'tes  sur  les  navires  et  les  azimuts.  Le  systeme  PolSAR  peut  detecter  des  navires  plus  petits 
que  les  systemes  a  polarisation  simple  ou  double. 

La  mission  de  la  Constellation  RADARSAT  (MCR)  comprendra  un  mode  CP  operationnel,  ce 
qui  devrait  etre  avantageux  pour  le  secteur  de  la  detection  des  navires.  II  est  recommande  de 
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considerer  la  possibilite  d’utiliser  le  mode  CP  pour  la  surveillance  a  grande  etendue,  en  particulier 
pour  la  detection  des  navires. 
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1  Introduction 


RADARSAT-2,  a  Canadian  satellite,  has  been  successfully  operating  in  space  since  late  2007.  It  carries 
a  synthetic  aperture  radar  (SAR)  that  provides  many  operating  modes,  including  the  Fine  Quad  (FQ) 
Mode,  which  provides  polarimetric  phase  and  amplitude  information. 

Improvements  in  ship  detection  capabilities  by  use  of  a  polarimetric  SAR  (PolSAR)  system  have 
previously  been  demonstrated  using  airborne  PolSAR  data  [1].  This  report  shows  that  RADARSAT-2 
FQ  imagery  also  provides  an  improvement  in  detection  performance.  It  was  shown  in  [1]  that  quad 
polarization  gives  the  best  performance,  followed  by  dual  polarization  systems,  which  in  turn  are  better 
than  single  polarization  systems. 

One  objective  of  this  report  is  to  describe  the  use  of  FQ  data  to  evaluate  the  performance  of  several 
types  of  SAR  systems:  single  polarization  HH,  dual  polarization,  horizontal  polarization  on  transmit 
and  horizontal  and  vertical  polarizations  on  receive  (HH  +  HV)  amplitude  only,  quad  polarization,  and 
compact  polarimetric  (CP).  The  decision  variable  used  for  the  detectors  is  based  on  the  likelihood  ratio 
test  with  the  Neyman-Pearson  criterion  as  defined  in  [1],  The  detection  performance  is  estimated  on  the 
basis  of  the  receiver  operating  characteristic  (ROC)  [2]. 

CP  radar  is  a  two  channel  approach  to  polarimetric  synthetic  aperture  radar  [4].  It  can  provide  a  wider 
swath  than  quad  polarization  because  it  uses  two  channels  rather  than  four.  For  maritime  surveillance 
applications,  wide  swath  coverage  is  essential  to  maximize  the  area  coverage.  For  example, 
RADARSAT  ScanSAR  Narrow  B  (SCNB)  Mode  images  acquired  using  dual  polarization  are  very 
favourable  for  ship  detection  [5]  because  of  the  wide  swath  width.  There  are  several  CP  configurations 
including  the  ti/4  mode  and  two  alternative  compact  modes.  In  the  tc/4  mode,  a  linearly  polarized  wave 
at  a  45°  orientation  is  transmitted  and  linear  H  and  V  polarized  waves  are  received  [6] [7].  In  the 
alternative  compact  modes,  a  circularly  polarized  wave  is  transmit  and  like  and  opposite  circularly 
(Compact  Circular  -  C,C),  or  linearly  H  and  V  polarized  waves  (Compact  Linear  -  CL)  are  received. 
The  CL  mode  is  also  referred  to  as  hybrid-polarity  [4]  [7]. 

In  this  study,  right-circular  transmit  and  linear  receive  (i.e.,  hybrid-polarity,  CL)  data  were  simulated 
using  RADARSAT-2  FQ  data.  The  application  of  PolSAR  ship  detection  algorithms  developed  using 
airborne  data  to  the  simulated  CL  data  was  investigated. 

SAR  data  of  the  Strait  of  Gibraltar,  which  contain  many  target  ships,  acquired  using  RADARSAT-2 
FQ,  were  used,  along  with  simultaneously  acquired  shore-based  automatic  identification  system  (AIS) 
data.  Detected  targets  were  validated  using  the  AIS  data. 

The  second  objective  is  to  optimize  the  detection  algorithms  based  on  the  likelihood  ratio  test. 

This  study  demonstrates  that  the  probability  density  function  (PDF)  of  the  ocean  decision  variable 
follows  approximately  a  K-distribution.  Thus,  the  algorithms  based  on  the  likelihood  ratio  test 
have  been  merged  with  a  K-distribution  constant  false  alarm  rate  (KCFAR)  detector.  Taking  an 
advantage  of  existing  KCFAR  detector  in  the  Analysts’  Detection  Support  System  (ADSS),  which 
was  developed  by  the  Defence  Science  and  Technology  Organization  (DSTO,  Australia)  [3],  the 
merged  algorithms  have  been  implemented  by  the  Defence  R&D  Canada  (DRDC). 
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In  this  report,  Section  2  provides  a  summary  of  RADARSAT-2’s  capabilities;  Section  3  presents  the 
methodologies  used  in  this  study;  Section  4  describes  CP  SAR;  Section  5  contains  a  discussion  of  the 
data  analysis  tools;  Section  6  summerizes  the  data  used  in  this  study;  Section  7  presents  the  results; 
Section  8  describes  the  KCFAR  Quad  detector  using  the  ADSS  and  the  initial  detection  results;  and 
Section  9  presents  the  conclusions. 


2 


DRDC  Ottawa  TM  2009-285 


2  RADARSAT-2  capabilities 


RADARS AT-2  is  a  C-band  (5.405  GHz)  SAR  satellite.  The  system  includes  polarimetric  SAR 
capability  and  supports  right-  and  left-look  imaging.  A  number  of  imaging  modes  provide  a  trade-off 
between  spatial  resolution  and  swath  coverage.  The  high-resolution,  narrow-swath  width  modes,  such 
as  Ultrafine,  have  a  nominal  resolution  of  3  m  with  a  swath  coverage  of  20  km.  In  contrast,  the  low- 
resolution,  wide-swath  width  modes,  such  as  ScanSAR  Wide,  have  a  nominal  resolution  of  100  m  with 
a  swath  coverage  of  500  km.  Certain  modes  are  available  in  a  choice  of  single,  dual  or  quad 
polarization.  The  operational  modes  and  the  coverage  are  illustrated  in  Figure  1  [8].  RADARSAT-2 
was  launched  in  December  2007;  it  is  owned  and  operated  by  MacDonald  Dettwiler  and  Associates  Ltd 

[9]- 

In  this  report,  we  focus  on  the  FQ  Mode.  In  FQ,  the  SAR  transmits  horizontal  and  vertical  polarized 
electromagnetic  waves  on  alternate  pulses  and  receives  the  scattered  signal  in  both  polarizations 
simultaneously.  The  quad  polarization  data  consist  of  amplitude  and  phase  information  of  measured 
scattered  fields  from  a  target.  The  additional  phase  information  aids  in  target  classification. 

The  FQ  Mode  is  intended  for  applications  that  require  polarimetric  information  (i.e.  amplitude  and 
phase)  of  a  target,  but  this  is  done  with  a  trade-off  of  narrow  swath  coverage.  The  FQ  Mode  can  cover 
an  incidence  angle  range  from  20°  to  41°  with  many  beam  positions;  each  beam  position  “covers” 
approximately  2°.  Only  one  beam  position  can  be  used  during  each  image  acquisition.  The  incidence 
angle  is  defined  as  the  angle  between  the  line  of  sight  from  the  radar  and  the  local  vertical  direction. 
The  nominal  resolution  is  5.4  m  x  7.9  m  (slant  range  x  azimuth)  for  single  look  complex  (SLC)  data 
and  the  scene  size  is  approximately  25  km  x  25  km. 
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Figure  1:  RADARSAT-2  operational  modes  and  coverage, 
(used  with  permission  [8]) 
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3  Overview  of  methodology 


3.1  Approach 

The  overall  approach  to  PolSAR  ship  detection  in  this  study  consists  of  many  steps  as  illustrated  in 
Figure  2.  The  most  important  step  is  the  ship  detection.  This  is  followed  by  verification  using  A1S  data. 
If  a  detected  candidate  target  has  corresponding  AIS  data,  then  the  candidate  target  is  confirmed  as  a 
ship.  If  the  candidate  target  does  not  have  AIS  data,  ambiguity  discrimination  using  phase  analysis  is 
used.  In  this  report,  we  focus  on  algorithms  developed  at  DRDC,  including  target  detection  and  azimuth 
ambiguity  analysis. 


Figure  2:  Overall  detection  approach. 
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3.2  LRT  ship  detection 


The  methodology  of  likelihood  ratio  test  (LRT)  ship  detection  using  polarimetric  SAR  data  has  been 
described  in  [1],  A  polarimetric  SAR  system  provides  observations  of  the  scattering  matrix  X(/,  j)  for 
each  pixel  (/,/)  in  an  image.  The  matrix  components  Snu(i,  j),  Sm{i,  j),  S  VH  (i,  j)  and  S  w(Uj)  give  the 
transformation  between  the  transmitted  and  scattered  electromagnetic  fields.  These  consist  of  the 
amplitude  and  phase  of  each  of  the  four  channels  of  the  polarimetric  data.  The  components  of  X(z,  j) 
can  be  written  as  a  vector: 

x(i,j)  =  [SHH(i,j)  SHV(i,j )  SVH(i,j )  Svv(i,j)]T  (1) 

where  the  superscript  T  is  the  transpose  operator. 

Ship  detection  is  a  binary  decision  problem.  In  the  fundamental  algorithms  of  polarimetric  SAR  ship 
detection,  statistical  decision  theory  is  applied  directly  to  the  components  of  the  scattering  matrix  to 
obtain  a  decision  variable.  A  likelihood  ratio  test  with  the  Neyman-Pearson  criterion  is  used  to  define  a 
pixel-based  detection  criterion.  Gaussian  distributions  for  each  of  the  scattering  matrix  components  are 
assumed  to  derive  the  formulation  while  measured  data  are  used  to  calculate  the  detection  variables. 

The  decision  variable  is  approximately  given  by  [1]: 

>77  for  a  ship 

x"(c;1)x  =  <  (2) 

<  77  for  ocean 

where  Ca  =  E(XXH)  is  the  covariance  matrix  of  the  ocean,  E(-)  is  the  expectation  operation,  H  is  the 
Hermitian  transpose  operator  and  77  is  the  designed  detection  threshold. 

As  shown  in  [1],  this  methodology  can  be  applied  to  dual  polarization  and  single  polarization  systems. 
In  a  dual  polarization  system,  the  vector  X(zj)  has  two  elements  and  in  a  single  polarization  system, 
X(zj)  has  one  element. 

A  detection  threshold  is  calculated  based  on  a  specified  false  alarm  rate  (FAR)  which  is  interpreted  as 
the  probability  of  false  alarm  (. PFA)  by  using  a  combination  of  extrapolation  and  curve  fitting  methods 
[11]. 

As  shown  in  [1],  a  statistical  decision  variable  is  calculated  for  each  pixel  of  the  selected  ocean  areas  in 
the  image  using  Equation  (2).  This  is  followed  by  the  calculation  of  the  median  value.  The  median  is 
used  as  an  initial  threshold  value  because  it  is  less  affected  than  the  mean  by  the  presence  of  large  pixel 
amplitudes  if  ships  are  present.  Finally,  the  other  thresholds  are  calculated  in  increments  of  half  of  the 
median.  A  curve  of  PFA  as  a  function  of  these  thresholds  is  then  determined.  The  number  of  pixels  is 
limited  in  the  selected  ocean  area,  consequently,  the  PFA  obtained  from  the  measured  data  is  limited.  In 
order  to  obtain  a  specified  PFA,  such  as  10’6  or  10"8,  an  extrapolation  method  is  used.  A  2nd  order 
polynomial  is  fitted  to  the  data  points  (PFA  versus  threshold)  using  a  minimum  mean  square  error 
method.  The  2nd  order  polynomial  equation  is  then  used  to  calculate  the  threshold  for  the  desired  PFA 
[11]. 
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Estimation  of  detection  performance  is  provided  by  ROC  curves,  which  are  a  plot  of  the  probability  of 
missed  detection  ( Pmd )  versus  probability  of  false  alarm  ( Pfa ).  These  probabilities  are  determined  by 
varying  the  detection  threshold. 

The  algorithms  described  here  are  for  pixel-based  detection;  complete  ship  detection  is  performed  after 
pixel  detection  by  using  a  clustering  technique.  (Examples  can  be  found  in  [3].)  The  PFA  specified  in 
this  report  is  for  pixel  detection.  The  application  of  clustering  will  reduce  the  overall  PFA  for  ship 
detection  because  more  than  one  detected  pixel  is  required  to  declare  the  presence  of  a  target. 


3.3  Azimuth  ambiguity  analysis 

In  a  radar  image,  azimuth  ambiguities  arise  due  to  discrete  sampling  of  the  Doppler  spectrum  at  the 
interval  of  the  pulse  repetition  frequency  (PRF  or  fp)  [12],  Since  the  Doppler  spectrum  repeats  at  fp 
intervals,  the  sampled  signal  components  outside  of  this  frequency  interval  will  be  aliased  into  the 
central  region  of  the  spectrum.  Ambiguities  are  often  detected  as  candidate  targets  in  the  maritime 
environment  since  the  clutter  (ocean)  signal  may  be  low.  However,  azimuth  ambiguities  can  be 
distinguished  from  the  main  target  image  by  using  analysis  of  the  phase. 

Under  the  reciprocity  assumption,  the  HV  and  VH  channels  are  identical  (i.e.,  SHv  =  SVH  for  both 
amplitude  and  phase).  Also,  the  phase  difference  between  samples  in  the  HV  and  VH  channels  is  zero 
for  targets  and  even-order  azimuth  ambiguities,  but  is  n  for  odd-order  azimuth  ambiguities  [13].  This 
phase  difference  between  the  HV  and  VH  channels  can  be  used  to  discriminate  between  targets  and 
odd-order  azimuth  ambiguities. 

The  azimuth  shift  of  the  ambiguity  location  relative  to  the  true  target  location  can  be  determined  [14]: 


AX 


amb 


2VS  { 


—  COS  {(p) 
0)  ) 


(3) 


where  X  is  the  radar  wavelength,  R  is  the  slant  range  at  which  the  target  is  broadside  to  the  radar,  Vs  is 
the  spacecraft  velocity,  coe  is  the  earth  rotation  speed  (2rc/day),  eo  is  the  spacecraft  orbital  rate  (343 
orbits  in  24  days)  and  cp  is  the  inclination  of  the  orbit  of  the  spacecraft. 
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4  Simulated  compact  polarimetric  SAR 


CP  is  a  type  of  dual  channel  synthetic  aperture  radar  in  which  one  polarization  is  transmitted,  but  two 
orthogonal  polarizations  are  received,  providing  both  phase  and  amplitude  information  of  the  measured 
backscattered  fields.  Target  backscatter  information  is  available  through  access  to  Stokes  parameters. 
However,  a  CP  system  cannot  provide  “pure”  HH  or  VY  backscatter  information  and  cannot  be 
considered  to  be  a  substitute  for  a  quad-pol  system. 

A  quad-pol  system  is  expected  to  provide  the  best  ship  detection  results  compared  to  other  SAR 
systems  at  the  same  resolution  [1]  since  it  contains  full  characterization  of  the  backscattered  field,  but  at 
the  expense  of  a  narrow  swath.  A  CP  system  is  expected  to  benefit  ship  detection  when  compared  to 
single-polarization  and  conventional  dual-polarization  systems  since  it  contains  relative  phase 
information  between  the  two  receive  channels.  The  CP  system  requires  only  one-half  of  the  data 
volume  per  pixel  compared  to  a  quad-pol  system,  so  it  should  be  possible  to  achieve  a  wider  swath 
width  than  a  quad-pol  system  (i.e.,  double  the  swath  width).  The  trade-off  between  the  information, 
swath  width  and  system  complexity  for  the  same  resolution  is  clearly  demonstrated  in  the  polarization 
information  hierarchy  table  (Table  1). 


Table  1:  Polarization  information  hierarchy. 


P 

o 

H 

fD 

ffl 

Ui 

s 

I 

09 

fD 


81 


el 

0 

p 


Mode 

Transmit 

Observe 

Swath 

SINGLE 

POLARIZATION 

H  or  V 

SHHov  Synov 

S hv  0^  $ w 

Wide 

ALTERNATING 

POLARIZATION 

H  and  V 

(Shh,  Syy) 

(no  phase) 

Wide 

DUAL 

POLARIZATION 

ON  RECEIVE 

H  or  V 

(Shh,  Shv)  or 
(■ Sw j  Syn) 

Wide 

COMPACT 
POLARIMETRY 
(CL  MODE) 

Circular 

(Shh~  j*SW)/21/2 
(SvH~  jSyy)/2U2 

Wide 

QUAD 

POLARIZATION 

H  and  V 

Shh,  Sm’, 

SvH ,  &  Syy 

Narrow 

The  various  possible  CP  modes  have  been  presented  in  Section  1.  Comparison  of  these  CP  SAR  modes 
has  been  discussed  in  the  literature  e.g.  [4][7][15].  For  example,  it  has  been  shown  that  the  well-known 
Cloude-Pottier  entropy-alpha  classification  can  be  performed  using  the  CL  system  and  the  n/4  mode  is 
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almost  as  good  as  full  quad-pol  for  crop  classification  [15].  Also,  circular  polarization  on  transmission 
is  preferred  when  there  is  Faraday  rotation  in  the  propagation  medium  [7].  In  addition,  the  architecture 
of  a  CL  system  is  relatively  simple  [4], 

In  the  context  of  this  study,  the  received  components  of  a  CL  system  consist  mostly  of  Shh  and  Syr 
because  cross-pol  components  are  usually  small.  It  has  been  shown  that  a  dual-pol  system  that  uses 
HH  and  VV  polarization  combinations  provides  good  ship  detection  performance  [1],  It  has  also  been 
shown  that  the  detection  performance  of  HH  and  VV  is  usually  better  than  for  conventional  dual-pol 
systems  with  single  transmit  polarization  (HH  and  HV,  or  VV  and  VH).  Therefore,  it  is  expected  that  a 
CL  system  should  also  be  capable  of  better  ship  detection  performance  than  conventional  dual-pol 
systems  and  this  is  possible  with  a  wider  swath  than  a  quad-pol  system. 


Right-circular  transmit  and  linear  receive  SAR  can  be  simulated  from  RADARSAT-2  FQ  data  as 
shown  in  [4]  and  [16]: 


e  e 

^  HH  ^  HV 

1 

ll 

1 

~sHH 

jS  HV 

C 

^  HV  ^VV 

_-j\ 

'77 

SVH 

-  jSvv  _ 

where  J  is  the  Jones  vector.  The  scattering  matrix  components  are  taken  from  the  measured 
RADARSAT-2  FQ  data. 

If  we  also  consider  system  noise  effects,  then  the  terms  of  the  X  vector  for  the  CL  system  become: 


X 


(^+Ar1)-7CV+V2) 

.(SpH+N3)-jXSrv+Ni)_ 


(5) 


where  Ni,  N?,  N3  and  N4  are  the  receiver  noise  in  each  channel  [16].  The  signal  and  noise  components 
in  each  channel  are  assumed  to  be  Gaussian  and  independent,  with  identical  noise  power  E([/V]2).  Using 
the  independence  properties  of  the  signal  and  noise  components,  it  can  be  shown  that  the  system  power 
in  each  channel  is: 

E(  |X,|!  )  =  |-  E(  |SHH  p  +  j  S  HH  ■  S’HV  -  j  S  HV  ■  Sh„  +  |SHV  |2  +  2- 1 N  |2 )  W 

E(  |X2|2  )  =  2  e(  |Svh[  +  jSVH  -Stv  -  jS„  -S™  +|SW|!  +2-|  N  |2)  (7) 


where  E  is  the  expectation  operator.  An  actual  CP  system  would  have  half  the  noise  power  in  each 
channel  compared  to  a  CP  system  simulated  using  RADARSAT-2  FQ  data.  However,  the  main  terms 
contributing  in  each  channel  are  co-pol  terms,  i.e.,  SHH  and  Svv,  which  are  usually  strong  signals. 

The  system  performance  is  approximately  determined  by  the  effective  target  to  clutter  ratio  of  the 
signal  components.  The  effect  of  the  noise  is  to  increase  the  apparent  power  of  each  signal  component. 


DRDC  Ottawa  TM  2009-285 


Typical  system  noise  data  for  the  image  which  was  used  is  illustrated  in  Figure  3,  which  shows  the 
levels  of  co-pol  and  cross-pol  clutter  relative  to  the  system  noise  floor.  The  noise,  taken  from  the 
product  meta  data,  is  expressed  as  noise  equivalent  sigma  zero  (NESZ).  The  NESZ  graph  of  Figure  3 
represents  the  noise  in  the  ETH  channel  whose  power  is  approximately  the  same  as  that  in  other 
channels.  In  the  top  graph,  the  backscatter  (sigmaO,  or  a0)  of  each  channel  (HH,  HV,  VH,  and  VV)  is 
shown  for  the  area  of  interest,  while  in  the  bottom  graph,  the  o°  is  shown  across  the  entire  range  for  the 
same  azimuth  area  as  the  top  graph. 


RSAT2,2C0e-07-29,HH  VV  HV  VH,Q4,Asc,Ls1art :  1467,Lstop  :  26S5,Pstart :  364,Pstop  :  1492 


HH-Pol  -  VV-Pol  HV-Pol  VH-Pol  NESZ 


RSAT2.2008  07  29.HH  VV  HV  VH,Q4,Asc,L3lart :  1467,Lslop  :  2655,Full  Swath 


Figure  3:  Graphs  of  clutter  and  noise  levels  for  R-2  FQ  data. 


In  the  case  of  Figure  3,  the  clutter  to  noise  ratio  for  the  co-pol  components  of  the  ocean  clutter  is  about 
30  dB,  while  the  clutter  for  the  cross-pol  components  is  only  about  6  dB  above  the  noise  floor.  The 
noise  power  in  each  channel  of  the  simulated  CL  product  is  increased  by  3  dB  compared  to  single 
channel  because  noise  from  both  co-pol  and  cross-pol  channels  is  present.  The  clutter  to  noise  ratio  at 
each  linear  polarization  component  will  thus  be  decreased  by  3  dB. 

Figure  3  also  shows  a  number  of  targets  in  the  cross-pol  transect  that  have  target  to  clutter  ratios 
ranging  from  7  to  16  dB.  The  cross-pol  targets,  which  are  13  to  22  dB  above  the  noise,  will  have  their 
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apparent  levels  raised  by  0.2  to  0.03  dB  respectively.  The  co-pol  clutter  apparent  level  will  be 
increased  by  only  about  0.004  dB  compared  to  the  level  without  noise  because  it  is  30  dB  above  the 
noise.  Thus,  the  target  to  clutter  ratio  for  the  cross-pol  components  will  be  lowered  by  less  than  1  dB 
compared  to  the  value  without  noise  while  the  target  to  clutter  ratio  of  the  co-pol  components  will  be 
degraded  by  a  negligible  amount. 

For  the  simulated  CL  system,  Equations  (6)  and  (7)  show  that  each  channel  will  contain  a  co-pol 
component  as  well  as  cross-pol  components  and  a  doubled  noise  power.  The  dominant  component  of 
each  channel  is  the  co-pol  component.  The  co-pol  clutter-to-noise  ratio  will  now  be  about  27  dB 
because  of  the  doubled  noise  and  the  apparent  level  of  the  co-pol  component  will  be  increased  by  about 
0.009  dB,  which  is  still  negligible.  The  levels  of  target  signals  will  also  be  unaffected  by  the  noise 
because  they  are  much  larger  than  the  clutter.  Thus,  the  effective  target-to-clutter  ratio  of  each  channel 
of  the  simulated  CL  system  will  be  negligibly  affected  by  the  presence  of  the  noise. 

It  is  concluded  that  using  R-2  FQ  Mode  data  to  simulate  the  CL  system  will  provide  a  realistic 
performance  estimate  which  is  unaffected  by  the  system  noise  because  the  received  signal  channels  are 
dominant  by  the  co-pol  components. 
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5  Data  analysis  tools 


5.1  PolSAR  tools 

Most  of  the  data  analysis  tools  used  in  this  study  were  developed  at  DRDC  Ottawa.  These  tools  include 
the  Likelihood  Ratio  Test  PolSAR  Ship  Detector,  polarization  synthesis,  phase  analysis  and  target 
signature  analysis,  as  described  in  Section  3  and  Figure  2.  The  tools  were  developed  in  Matlab™. 


5.2  IA  Pro 

In  this  report,  some  of  the  results  are  displayed  using  Image  Analyst  Pro  (IA  Pro)  which  was  developed 
at  DRDC  Ottawa  for  validation  and  demonstration  of  new  algorithms  for  GEOINT  [17].  The  focus  of 
IA  Pro  development  is  on  tools  that  introduce  automation,  feature  extraction  and  image 
combination/fusion  to  the  exploitation  of  SAR  and  EO  imagery.  IA  Pro  provides  support  for  single-  and 
multi-channel  SAR  and  EO  imagery,  and  permits  the  image  analyst  to  incoiporate  geospatial 
information  (i.e.,  thematic  vector  layers  and  digital  maps)  for  contextual  awareness. 

IA  Pro  was  built  using  C++,  Python,  OpenEV  and  the  Geospatial  Data  Abstraction  Library  (GDAL). 
Python  is  an  open-source  object-oriented  interpreted  programming  language,  OpenEV  is  an  open- 
source  library  of  raster  and  vector  classes  and  functions,  and  GDAL  is  a  translator  library  for  raster 
geospatial  data  formats.  IA  Pro  provides  a  multitude  of  tools  for  the  manipulation  and  exploitation  of 
multi-sensor  and  multi-temporal  imagery  for  GEOINT.  A  detailed  description  of  IA  Pro  can  be  found  in 
[17]. 


5.3  ADSS 

The  Analysts’  Detection  Support  System  is  an  extensible  suite  of  image  processing  algorithms  for 
GEOINT  applications  -  including  automatic  target  detection  (ATD)  for  broad-area  search  for  small 
targets  -  with  SAR  and  EO/IR  sensors  [3].  It  was  developed  by  DSTO,  with  additional  contributions 
from  the  Defence  Science  and  Technology  Laboratory  (Dstl)  and  DRDC.  ADSS  includes  algorithms  for 
ATD,  change  detection,  and  region  classification  for  both  maritime  and  land  areas,  in  addition  to  a 
number  of  support  and  utility  modules.  The  flexibility  of  ADSS  arises  from  a  system  architecture  that 
supports  the  combination  of  algorithms  in  a  user-specified  structure.  (The  numerous  algorithms  and 
utilities  that  are  provided  by  ADSS  are  implemented  as  modules  that  can  be  combined  by  using 
Command  and  Data  Language  (CDL)  scripts.)  The  power  of  ADSS  originates  in  a  framework  that 
provides  for  synchronization,  control,  and  data  flow  between  the  various  modules. 

In  the  second  part  of  this  study,  a  popular  and  widely  used  ATD  algorithm  the  K-distributed  constant 
false  alarm  rate  (K-CFAR)  algorithm  (KCFAR  module)  implemented  in  ADSS  [18]  [19]  was  applied  to 
the  data.  The  KCFAR  module  uses  the  K-distribution  to  model  the  probability  density  function  (PDF) 
of  inhomogeneous  ocean  clutter  and  system  noise  in  high-resolution  SAR  images.  The  KCFAR  module 
locally  estimates  the  parameters  for  the  K-distribution  and  assigns  detections  to  bright  points  appearing 
in  the  tail  of  the  distribution. 
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6  Data 


6.1  Polarimetric  data 

The  polarimetric  data  used  in  this  study  were  acquired  in  RADARSAT-2  FQ  Mode  on  29  July 
2008  (PDS  00278000).  The  Strait  of  Gibraltar  was  chosen  since  there  are  many  ships  in  this 
region  and  there  are  validation  data  available  from  AIS.  The  data  were  acquired  on  an  ascending 
pass  and  the  nominal  incidence  angle  range  is  from  22°  to  24°.  The  signal  data  were  processed  to 
a  calibrated  image  by  MDA.  Radar  backscatter  coefficient  (a0)  values  were  used  in  this  analysis. 


6.2  AIS  data 

AIS  is  an  internationally  mandated  ship-based  transponder  system  [19]  that  is  used  by  vessels  to 
inform  one  another  about  their  activity.  Information  available  includes  identity,  location,  speed, 
heading,  type.  International  Maritime  Organization  (IMO)  number,  and  Maritime  Mobile  Service 
Identity  (MMSI)  number.  This  system  is  designed  primarily  for  safety  of  life  at  sea  (SOLAS), 
collection  avoidance  and  traffic  management.  AIS  operates  in  the  VHF  maritime  band  using  a 
Self-Organized  Time  Domain  Multiplex  (SOTDMA)  protocol  that  coordinates  communication 
between  vessels  within  line-of-sight;  shore-based,  air-based,  or  space-based  receivers  can  also 
receive  this  information.  Although  this  cooperative  identification  system  is  intended  for  collision 
avoidance,  it  is  available  to  assist  maritime  wide  area  surveillance  by  providing  vessel  identities. 

AIS  data  were  acquired  along  with  the  RADARSAT-2  data  and  are  used  to  verify  identities  and 
locations  of  ships  detected  in  this  study.  The  AIS  data  were  obtained  from  the  Maritime  Safety 
and  Security  Information  System  (MSSIS)  which  was  initiated  by  the  US  Coast  Guard. 
International  contributors  gain  access  to  the  database  which  is  generated  from  real  time  reports 
from  shore-based  receivers.  The  MSSIS  data  used  for  this  report  were  provided  by  DRDC 
Atlantic  and  processed  at  DRDC  Ottawa  [5]. 

Ten  ships  for  which  AIS  data  are  available  have  been  used  in  this  and  other  separate  studies. 
Eight  of  them  have  been  selected  for  detection  performance  analysis  (Section  7.5).  Ship 
information,  such  as  length,  width,  speed  and  type  can  be  extracted  from  the  AIS  data,  as  listed  in 
Table  2. 
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Table  2:  Extracted  ship  information  from  AIS  data  (July  29,  2008). 


Ship 

Speed 

(knots) 

Width 

(m) 

Length 

(m) 

Draught 

(m) 

Type 

1 

0 

30 

258 

7.8 

Cargo 

2 

1.1 

42 

332 

11.8 

Tanker 

3 

0.3 

32 

181 

10.1 

Tanker 

4 

0.3 

49 

288 

10.5 

Tanker 

5 

1.1 

6 

18 

4.2 

Tug 

6 

16.1 

30 

198 

10.6 

Cargo 

7 

0.7 

44 

277 

10.2 

Tanker 

8 

15.2 

60 

333 

11.5 

Tanker 

9 

10.35 

25 

184 

7 

Cargo 

10 

0.22 

26 

198 

8 

Tanker 
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7  Results 


7.1  Quad  polarization  image 

A  quad  polarization  image  can  be  displayed  in  many  ways,  for  example  as  a  single  channel  in 
grey  scale  or  a  linear  combination  of  four  channels  in  an  RGB  color  image.  In  this  study,  the 
Pauli-basis  RGB  display  is  used  since  it  permits  the  extraction  of  physical  information.  In  a 
Pauli-basis  image,  the  Shh  +  $vv  component  represents  single  (odd)  bounce  scattering,  while  Shh~ 
Syy  represents  double  (even)  bounce,  and  Sm-  +  Svh  represents  volume  scattering.  The  physical 
interpretation  of  these  three  basic  scattering  mechanisms  is  illustrated  in  Figure  4. 


SINGLE  BOUNCE  DOUBLE  BOUNCE 

SCATTERING  SCATTERING 


VOLUME 


HV  +  VH 


Figure  4:  Physical  interpretation  of  three  basic  scattering  mechanisms. 


The  Strait  of  Gibraltar  image  which  is  used  as  a  test  image  is  displayed  in  the  Pauli-basis  in 
Figure  5.  In  this  image,  the  magnitude  of  ax|HH  -  VV|  is  displayed  in  red,  the  magnitude  of 
bx|HV  +  VH|  is  displayed  in  green,  and  the  magnitude  of  cx|HH  +  VV|  is  displayed  in  blue.  The 
coefficients  a,  b,  and  c  are  selected  to  be  2,  3  and  1  respectively.  The  ocean  tends  to  exhibit  single 
bounce  scattering,  while  ships  tend  to  have  a  significant  amount  of  double  bounce  scattering. 
Vegetation  tends  to  have  significant  volume  scattering  in  addition  to  surface  scattering. 
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Figure  5:  Pauli-basis  RGB  image  of  the  Strait  of  Gibraltar  in  a  geodetic  projection. 

Original  RADARS AT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates 
Ltd.,  2008  -  All  Rights  Reserved. 
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7.2  LRT  ship  detection 

The  Likelihood  Ratio  Test  (LRT)  Ship  Detection  software  was  applied  to  the  test  image  and  the 
detection  threshold  was  calculated  based  on  a  specified  probability  of  false  alarm  (Pfa)  using  a 
combination  of  the  extrapolation  and  curve  fitting  methods,  as  described  in  Section  3  and  in  [1 1]. 

In  order  to  calculate  the  threshold,  we  need  to  select  an  area  of  the  ocean.  In  this  example,  the 
ocean  area  (range  pixels:  1000  to  1300  and  azimuth  pixels:  1600  to  2600)  consists  of 
approximately  3xl06  pixels.  The  detection  threshold  is  approximately  16  dB  and  is  calculated 
with  the  PFA  set  to  10'6  for  detection  in  the  quad-pol  image.  The  images  of  the  ocean  decision 
variable  used  to  calculate  the  threshold  are  illustrated  in  Figure  6a)  and  b).  In  Figure  6a),  we  can 
see  that  there  is  a  candidate  target  in  the  selected  ocean  region. 

In  order  to  reduce  or  even  eliminate  the  effect  of  pixels  with  high  decision  variable  values  (i.e. 
potential  targets  or  land)  on  the  calculation  of  the  threshold,  two  steps  were  involved  in  the 
calculation.  The  detection  variable  is  estimated  from  the  selected  ocean.  Then,  the  threshold  was 
calculated  using  only  pixels  which  have  decision  variable  values  less  than  3  times  the  median  of 
the  decision  variable,  as  illustrated  in  Figure  6b),  which  indicates  that  candidate  target  was 
excluded  (appears  black  in  the  image  since  the  values  are  zero).  Figure  6b)  is  a  zoom  in  of  the 
target  location  in  Figure  6a). 

In  this  section,  the  focus  is  on  the  application  of  the  detection  algorithms.  Optimization  of  the 
detection  is  described  later  in  this  document. 

The  detection  results  are  shown  in  Figure  7  where  pixels  having  values  smaller  than  the  detection 
threshold  were  assigned  to  be  zero  and  the  pixel  values  were  retained  for  pixels  having  values 
greater  than  the  threshold.  In  Figure  7,  the  green  boxes  indicate  locations  where  ship  signatures 
are  predicted  to  be  based  on  the  MSSIS  AIS  data  and  the  image  acquisition  parameters.  Based  on 
the  available  AIS  data,  34  out  of  36  ships  were  detected  and  the  smallest  of  these  had  a  length  of 
18  m  and  a  speed  of  0.59  m/s.  Two  ships  identified  in  the  AIS  data  were  not  detected  by  the 
PolSAR  detector.  One  of  them  was  very  close  to  land,  as  indicated  in  this  figure. 

An  example  of  a  detected  unknown  ship  target  is  indicated  as  Tu.  The  azimuth  ambiguity  of  this 
target  was  also  detected  as  a  candidate  target  (Au).  It  has  been  found  that  azimuth  ambiguities  can 
be  distinguished  from  targets  because  of  their  phase  characteristics,  as  discussed  in  the  next 
section. 
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Figure  6:  a)  Image  of  ocean  decision  variable;  b)  zoom  in  of  an  excluded  target 
area  from  ocean  samples  in  a). 

Original  RADARS AT-2  Data  and  Products  ©MacDonald,  Dettwiler  and 
Associates  Ltd.,  2008  -  All  Rights  Reserved. 
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□  AIS  □  Target  of  interest  (see  Figure  12  and  13) 

Figure  7:  Quad-pol  detection  results  of  the  Strait  of  Gibraltar  in  a  geodetic  projection. 

Original  RADARS AT-2  data  and  products  ©  MacDONALD,  DETTWILER  AND  ASSOCIATES 
LTD.  (2008)  -  All  Rights  Reserved. 
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7.3  Phase  analysis 


The  correlations  of  and  phase  difference  between  channels  corresponding  to  different  polarization 
were  investigated.  The  correlation  and  phase  difference  between  channels  were  calculated  using  a 
moving  average  with  a  window  size  of  5x5  pixels.  (Different  window  sizes  have  not  been  tested 
in  the  present  study.)  The  absolute  values  of  correlation  coefficients  for  ocean  are  calculated 
using  the  selected  ocean  area  (range  pixels  1000  to  1300  and  azimuth  pixels  1900  to  2600  in 
Figure  6a) )  and  the  results  are  given  in  Table  3. 


Table  3:  Absolute  correlation  coefficients  for  ocean. 


HH 

HV 

VH 

VV 

HH 

1.00 

0.05 

0.07 

0.97 

HV 

0.05 

1.00 

0.46 

0.04 

VH 

0.07 

0.46 

1.00 

0.07 

VV 

0.97 

0.04 

0.07 

1.00 

The  correlation  of  and  phase  difference  between  the  HH  and  VV  channels  are  shown  in  Figure  8 
and  9  respectively.  In  the  phase  difference  image  (Figure  9),  the  phase  angle  is  indicated  in 
degrees  in  the  colorbar. 

From  the  results,  we  can  see  that  the  image  of  the  ocean  is  highly  correlated  between  HH  and  VV 
channels.  This  is  consistent  with  our  previous  studies  and  results  in  the  literature  [13][20].  The 
high  correlation  of  the  ocean  aids  in  observation  of  candidate  targets,  consequently  it  aids  in 
discriminating  the  target  and  the  “ghost”  images.  These  “ghost”  images  in  the  range  direction  are 
due  to  the  payload  performance  for  high  resolution  modes.  Since  these  data  were  acquired,  MDA 
has  implemented  modifications  to  significantly  reduce  these  “ghost”  images.  Two  examples  of 
“ghost”  images  are  indicated  in  Figure  8. 
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Figure  8:  Correlation  coefficient  of  HH  and  VV for  test  area. 
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phase  of  HH  and  W 


Figure  9:  Phase  difference  between  HFl  and  VV for  test  area. 


The  correlation  of  and  phase  difference  between  the  HV  and  VH  channels  are  shown  in  Figure  1 0 
and  11,  respectively.  These  figures  indicate  that  the  targets  are  highly  correlated,  but  the  azimuth 
ambiguities  have  a  phase  difference  of  180°,  as  expected  [13].  As  previously  indicated,  this  phase 
difference  between  the  HV  and  VH  channels  can  be  used  to  aid  in  distinguishing  between  targets 
and  azimuth  ambiguities. 

Figure  12  gives  examples  of  correlation  of  and  phase  difference  for  a  few  targets  whose  locations 
are  indicated  in  Figure  7. 
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Figure  10:  Correlation  coefficient  of  HV  and  VH  for  test  area. 
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phase  of  HV  and  VH 


Figure  11:  Phase  difference  between  HV  and  VH  for  test  area. 


Figure  12:  a)  Correlation  coefficient  of  HH  and  VV;  b)  correlation  coefficient  of  HV  and  VH;  c) 
phase  difference  between  HV  and  VH. 
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The  correlation  of  the  HH  and  VV  channels  as  shown  in  Figure  12a)  indicates  that  the  ocean  is 
highly  correlated,  but  targets  are  not.  On  the  other  hand,  the  correlation  of  HV  and  VH  channels 
shown  in  Figure  12b)  indicates  that  the  targets  are  highly  correlated,  and  the  phase  difference  is 
approximately  180°,  as  shown  in  Figure  12c).  The  phase  analysis  results  confirm  that  detected 
candidate  targets  are  not  processing  artefacts.  The  signatures  of  some  ships  can  be  seen. 

Figure  13a)  shows  the  image  of  the  correlation  of  the  HV  and  VH  channels  of  a  target,  Figure 
13b)  shows  the  phase  difference  between  the  HV  and  VH  channels  of  the  ambiguity  of  this  target, 
and  Figure  13c)  shows  the  ambiguity  of  the  target.  This  ambiguity  was  detected  as  a  candidate 
target.  However,  it  was  distinguished  from  the  actual  targets  by  the  phase  analysis. 


Figure  13:  Example  phase  analysis  results  ofHV  and  VH.  a)  correlation  coefficient  for  a  ship:  b) 
phase  difference  for  a  target  and  an  ambiguity;  c)  phase  difference  for  an  ambiguity. 


Azimuth  ambiguities  can  also  be  clearly  observed  in  the  |HV-VH|  image.  Based  on  the  phase 
properties  of  HV  and  VH  channels  described  in  Section  3.3,  when  the  |HV-VH|  image  is  formed, 
the  amplitudes  of  the  odd-order  azimuth  ambiguities  are  enhanced,  while  the  amplitudes  of  the 
targets  and  even-order  azimuth  ambiguities  are  “cancelled”. 

The  |HV-VH|  image  of  the  scene  is  shown  in  Figure  14a).  Some  azimuth  ambiguities  can  be 
seen.  Examples  of  ambiguities  of  regions  on  the  land  are  indicated.  In  the  |HV-VH|  image,  the 
impact  of  the  system  noise  floor  higher  than  the  signal  is  visible  in  the  very  low  backscatter 
regions.  To  aid  in  observation  of  ambiguities  in  the  |HV-VH|  image,  the  |HV|  image  is  also 
shown  in  Figure  14b). 

According  to  Equation  (3),  the  distance  between  the  ambiguity  Ay  and  the  target  Tu  is  estimated 
to  be  approximately  4.4  km.  Similarly,  the  azimuth  ambiguities  (Aiand)  corresponding  to  the 
structures  on  the  land  which  occur  at  smaller  slant  ranges  are  found  to  be  displaced  approximately 
4.3  km  from  the  actual  position  (TLand)-  (See  also  Figure  13  for  the  target  Ty  and  the  ambiguity 
Ay.)  The  measurements  from  the  images  agree  with  these  calculated  values. 
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The  results  demonstrate  that  the  quad-pol  system  provides  improved  ship  detection  and  the  phase 
information  between  HV  and  VH  can  be  used  to  distinguish  azimuth  ambiguities  and  targets, 
consequently  aiding  in  reducing  the  false  alarm  rate.  The  relative  detection  performance  of 
various  polarization  combinations  will  be  considered  in  Section  7.5. 
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Figure  14:  a)  \HV  -  VH\;  b)  \HV\  image  of  the  Strait  of  Gibraltar. 

Original  RADARS AT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates 
Ltd.,  2008  -  All  Rights  Reserved. 
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7.4  Compact  PolSAR  ship  detection 

The  compact  polarimetric  SAR  data  were  simulated  using  measured  RADARSAT-2  FQ  Mode 
data  as  described  in  Section  4. 

The  Likelihood  Ratio  Test  Ship  Detection  software  was  applied  to  the  simulated  CL  data.  The 
detection  threshold  calculation  procedures  and  the  ocean  sample  area  are  the  same  as  described  in 
Section  7.2.  The  detection  threshold  is  approximately  13.4  dB  calculated  with  PFA  set  to 
approximately  10'6. 

The  detection  results  are  shown  in  Figure  15.  The  results  are  similar  to  those  from  the  quad-pol 
detection  shown  in  Figure  7. 
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AIS 


Figure  15:  Detection  results  for  simulated  CL  data  of  the  Strait  of  Gibraltar  in  a  geodetic 
projection. 


7.5  Ship  detection  comparison 

The  ship  detection  performance  for  RADARSAT-2  FQ  Mode  data  and  the  simulated  CL  data  as 
well  as  dual-pol  F1FI+F1V  amplitude  only  and  single-pol  HH  have  been  examined  and  the  ROC 
curves  determined.  The  HH  polarization  is  used  as  a  baseline  since  this  is  the  transmit-receive 
configuration  for  RADARSAT-1.  The  detection  performance  analysis  was  carried  out  for  eight 
ships  in  the  test  scene.  (See  Table  2.) 
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As  examples,  the  detection  performance  curves  for  ships  3  and  6,  are  presented  in  Figure  16a)  and 
17a),  respectively,  where  b)  and  c)  show  the  corresponding  |HH|  and  simulated  |X1|  images. 
Results  for  the  other  vessels  are  presented  in  Annex  A. 
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Figure  16:  a)  Detection  performance  for  ship  3;  b)\HH\  image  and  c)  simulated  X/  image. 

Original  RADARSAT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 
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b) 
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Figure  17:  a)  Detection  performance  for  ship  6;  b)\HH\  image  and  c)  simulated  \Xj\  image. 

Original  RADARSAT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 


In  general,  a  lower  Pfa  results  in  a  higher  Pmd,  or  a  lower  Pmd  results  in  a  higher  Pfa.  A  system 
having  the  lowest  Pfa  and  lowest  Pmd  provides  the  best  detection,  i.e.  the  lower  the  curve,  the 
better  the  performance.  Thus,  all  of  these  ROC  curves  clearly  show  the  advantage  (i.e.,  decreased 
probability  of  missed  detection,  decreased  probability  of  false  alarm)  of  a  quad  polarization 
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system,  which  uses  all  of  the  available  amplitude  and  phase  information  (triangles).  Dual 
polarization  systems  with  amplitude  only  provide  much  better  performance  than  single  channel 
HH  polarization  (asterisks). 

For  five  of  eight  cases  studied,  the  simulated  CL  data  (square)  provide  better  performance  than 
HH  +  F1V  with  amplitude  only  (circles).  For  the  other  three  cases,  HH  +  F1V  with  amplitude  only 
provides  slightly  better  performance  than  the  simulated  CL  system. 

By  arbitrarily  selecting  a  probability  of  false  alarm  of  Pfa  =  1 0"4,  the  relative  improvement  in  ship 
detectability  may  be  quantified  by  comparing  the  probabilities  of  missed  detection  Pud-  The 
relative  improvement  in  the  detection  performance  of  all  cases  studied  is  summarized  in  Table  4 
by  taking  the  ratio  of  the  observed  Pmd  to  that  of  quad  polarization  PmdQ-  We  observed  that  an 
approximately  seven- fold  decrease  in  the  probability  of  missed  detection  was  achieved  by  moving 
from  a  single  polarization  to  a  polarimetric  radar  system,  and  a  four-fold  decrease  by  moving  to  a 
simulated  CL  system.  In  some  cases,  the  single  channel  system  cannot  provide  such  low 
probability  of  missed  detection.  (See  Figure  16a).)  These  results  show  that  a  quad  polarization 
system  provides  the  best  detection  for  the  cases  studied,  but  a  CL  system  will  provide  much  better 
performance  than  a  single  channel  radar  system. 

A  statistical  analysis  of  the  detection  performance  of  these  systems  was  performed.  The  mean 
value  and  the  standard  deviation  of  the  mean  of  each  system  were  calculated  on  the  observed  Pmd 
values.  The  mean  values  of  the  PMd  and  the  standard  deviation  of  the  mean  values  are  normalized 
to  those  of  the  HH  system. 


The  uncertainty  in  the  normalized  mean  of  the  Pmd  is  given  by: 


SC  = 


(1  0  > 

2 

f  a  ^ 

2 

—  x  5a 

+ 

— rX  8b 

lb  j 

A  J 

(8) 


where  a  is  the  mean  of  individual  system,  Sa  is  the  corresponding  standard  deviation,  b  is  the 
mean  of  HH  system  and  Sb  is  the  corresponding  standard  deviation.  These  results  are 
summarized  in  Table  5. 
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Table  4:  Values  of  Pmd  far  Pfa  =  1  O'4  far  ships  considered. 


Ship  ID  # 

System 

Pmd 

P MtfP MlP 

1 

HH 

— 

HH+HV 

2-OxlO'1 

2.00 

CL 

L5xl0_1 

1.50 

Quad-pol 

l.OxlO'1 

1.00 

3 

HH 

— 

— 

HH+HV 

2-lxlO'1 

1.50 

CL 

2.2xl0_1 

1.57 

Quad-pol 

L4xl0_1 

1.00 

4 

HH 

— 

— 

HH+HV 

2.5xl0_1 

1.25 

CL 

3-OxlO'1 

1.50 

Quad-pol 

2-OxlO'1 

1.00 

6 

HH 

2-OxlO'1 

6.90 

HH+HV 

l.lxlO'1 

3.79 

CL 

4.6xl0'2 

1.59 

Quad-pol 

2.9xl0'2 

1.00 

7 

HH 

4.3xl0_1 

6.18 

HH+HV 

2-OxlO'1 

2.99 

CL 

l.OxlO'1 

1.49 

Quad-pol 

6.7xl0'2 

1.00 

8 

HH 

— 

— 

HH+HV 

L5xl0_1 

1.36 

CL 

L9xl0_1 

1.73 

Quad-pol 

l.lxlO'1 

1.00 

9 

HH 

4-OxlO'1 

4.08 

HH+HV 

2.9xl0_1 

2.96 

CL 

L3xl0_1 

1.33 

Quad-pol 

9.8xl0'2 

1.00 

10 

HH 

— 

— 

HH+HV 

3.2xl0_1 

2.91 

CL 

L4xl0_1 

1.23 

Quad-pol 

l.lxlO'1 

1.00 
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Table  5:  Statistical  analysis  results  of  detection  performance  for  PFA  =  1(T4. 


Mean  (PMD) 

std  of  mean 

C  (Mean/Mean™ )  ±  8C 

HH 

0.34 

0.07 

1.00  ±0.30 

HH+HV 

0.22 

0.02 

0.63  ±  0.15 

CL 

0.16 

0.03 

0.46  ±  0.13 

Quad-pol 

0.11 

0.02 

0.31  ±  0.08 

The  normalized  values  of  the  mean  of  the  PMD  for  the  quad  polarization,  simulated  CL  and 
HH+HV  systems  demonstrate,  from  another  perspective,  that  a  quad  polarization  system  provides 
a  significant  decrease  of  the  probability  of  missed  detection  from  that  for  a  single  polarization, 
while  a  CL  system  will  achieve  much  lower  probability  of  missed  detection  than  a  single  channel 
system. 

A  compact  polarimetric  system  should  be  capable  of  wider  swath  than  conventional  polarimetry 
having  similar  resolution  since  it  only  stores  half  of  the  data  volume  per  pixel.  Such  a  system  is 
expected  to  provide  better  detection  performance  than  a  conventional  dual  polarization  system. 
However,  more  investigation  is  required  to  verify  this  expectation.  In  contrast,  a  CL  system  could 
provide  higher  resolution  than  the  quad  polarization  system  while  maintaining  similar  swath 
coverage. 
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8  KCFAR  ship  detection 


Previous  sections  described  the  detection  performance  of  the  quad-pol  ship  detection  algorithms 
based  on  the  likelihood  ratio  test,  developed  using  airborne  PolSAR  data  to  the  RADARS AT-2 
FQ  data.  Further  optimization  of  these  algorithms  has  been  carried  out.  The  optimization, 
implementation  of  optimized  algorithms  and  initial  detection  results  are  presented  in  this  section. 


8.1  Estimation  of  K-distribution  parameters  for  quad-pol  data 

In  order  to  improve  the  quad-pol  ship  detection  described  in  Sections  3.3,  the  probability  density 
function  (PDF)  of  the  ocean  decision  variable  has  been  studied. 

The  K-distribution  for  ocean  samples  of  a  SAR  image  can  be  described  as  [18][22]: 


P(l)  =  j 


LvI 


1 


r(X)T(  y) 


-K„ 


(9) 


where  p(I)  is  the  probability  of  image  intensity  /,  (I)  is  the  mean  intensity,  L  is  the  equivalent 
number  of  statistically  independent  looks  (ENL)  in  the  scene,  v  is  the  order  parameter  and 
describes  the  local  modulation  of  the  mean  due  to  geophysical  processes  such  as  winds,  waves, 
and  currents,  T  is  the  gamma  function,  and  KV~L  is  the  modified  Bessel  function  of  the  second 
kind  of  order  v  -  L.  L  and  v  define  the  shape  of  the  K-distribution,  which  in  turn  are  used  in  the 
threshold  calculation  for  pixel-based  ship  detection  in  a  CFAR  detector. 

In  the  quad-pol  ship  detector,  the  image  intensity  /  is  the  decision  variable  calculated  using 
Equation  (2). 

In  general,  ENL  is  a  SAR  processor  characteristic.  However,  for  the  decision  variable,  which  is 
derived  from  the  combination  of  four  channels  of  single-look  products,  the  number  of  looks  L  is 
unknown.  Also,  the  order  number  v  is  unknown.  The  approach  which  we  have  taken  to  estimate 
these  two  parameters  is  as  follows. 

First,  we  choose  a  value  for  L  then  estimate  v,  then  plot  the  K-distributed  PDF  for  each  pair  of  L 
and  v,  then  test  how  well  the  estimated  K-distribution  function  fits  the  measured  data.  Two 
methods  were  used;  one  is  the  root  mean  square  error  and  the  other  is  the  yj  goodness-of-fit  test. 
If  the  tests  from  both  methods  agree  well,  then  the  L  and  v  are  found  acceptable.  The  L  is  varied 
from  1  to  8.  For  each  L,  v  is  varied  from  1  to  200.  The  increment  of  L  is  0.05  between  L  =  3.8  to 
4.1,  and  1  for  others.  For  the  ocean  samples  which  we  used  for  the  test,  it  was  found  that  L  is  4 
and  v  is  40.95. 

The  histogram  of  the  K-distributed  ocean  decision  variable  and  the  continuous  PDF  according  to 
Equation  (9)  are  shown  in  Figure  18.  The  results  indicate  that  the  K-distribution  is  a  good  model 
for  the  ocean  decision  variable  from  the  combination  of  four  channels. 


34 


DRDC  Ottawa  TM  2009-285 


Therefore,  the  detection  algorithms  based  on  the  likelihood  ratio  test  have  been  optimized  and  the 
decision  variable  has  been  input  into  a  KCFAR  detector.  The  KCFAR  detector  used  was  one 
developed  by  Dstl  in  ADSS  [19]:  the  algorithms  for  deriving  the  decision  variable  for  quad 
polarization  data  was  developed  by  DRDC  Ottawa  and  it  produced  the  image  input  to  the  KCFAR 
detector. 


L=4,  nuU=40.95.  mu=3.95 


Figure  18:  Probability  density  function  of  the  ocean  decision  variable. 


8.2  Operation  of  the  KCFAR  detector 

After  implementation  of  the  detection  algorithms  using  ADSS,  tests  were  performed  on  the  test 
image. 

In  order  to  avoid  significant  variation  of  ocean  statistics  over  large  areas  in  an  image,  the  decision 
variable  is  calculated  on  a  small  area  (tile)  basis.  The  entire  scene  is  divided  into  small  tiles:  i.e. 
10x10  tiles  for  FQ  images  since  the  coverage  is  only  25  x  25  km. 

In  order  to  eliminate  the  effect  of  ocean  pixels  with  high  decision  variable  values,  the  ocean 
sample  median  of  |HH|2  and  the  median  of  |HV|2  were  calculated,  then  the  |HHj2  and  |HV|2  values 
of  each  ocean  pixel  are  compared  to  these  medians.  If  the  intensity  of  each  pixel  meets  the 
condition  of  (|HH|2  <  5xmedian  |HH|2  and  |HV|2  <  5xmedian  |HV|2),  then  the  pixel  will  be  used  to 
calculate  the  decision  variable. 

The  decision  variable  for  each  pixel  in  an  image  is  calculated  and  sent  to  the  KCFAR  detector 
(KCFAR  Quad)  as  the  required  parameter  to  perform  the  pixel  detection. 
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The  straightforward  processing  chain  is  to  land  mask  the  images,  run  the  KCFARQuad 
algorithms,  cluster  the  results  into  candidate  targets,  and  display  the  results  in  1A  Pro. 

The  KCFAR  algorithm  in  ADSS  sets  a  guard  region  and  an  annulus  in  estimating  the  K- 
distribution  parameters  of  the  background  without  contamination  from  potential  target  pixels. 
Several  sets  of  values  for  the  annulus  parameters  were  tried  to  determine  the  behaviour  of  the 
algorithm.  Those  values  were  chosen  to  maximize  the  detection  of  targets  declared  on  the  basis  of 
A1S  data  and  the  visual  evaluation. 

For  this  study,  Pfa  =  10"10  is  used.  Other  parameters,  including  multi- look,  correlation  and  clutter 
estimation  method,  are  set  to  default  values. 


8.3  Initial  detection  results 


Initial  detection  results  from  the  KCFAR  Quad  detector  are  illustrated  in  Figure  19.  The  image 
(acquired  on  29  July  2008)  is  the  same  as  used  in  earlier  studies  described  in  this  report.  This 
image  shows  the  decision  variable  image  with  a  land  mask  based  on  the  World  Vector  Shoreline 
(WVS)  applied.  Detection  results  are  indicated  as  red  squares  while  the  A1S  results  are  green 
squares.  Thirty-six  ships  were  declared  by  the  AIS.  Two-hundred-forty  candidate  targets  were 
declared  by  the  KCFAR  Quad,  including  30  ships  declared  by  the  AIS.  The  smallest  ship 
declared  by  both  AIS  and  the  KCFAR  Quad  detectors  had  a  length  of  18  m.  It  is  noted  that  six  of 
the  ships  declared  by  the  AIS  appear  to  have  been  eliminated  by  application  of  the  land  mask. 

Some  azimuth  ambiguities  were  detected  as  candidate  targets.  The  azimuth  ambiguities  can  be 
distinguished  from  the  targets  according  to  Equation  (3)  in  section  3.3.  Examples  of  azimuth 
ambiguities  which  were  detected  as  candidate  targets  are  shown  in  Figure  19  (the  decision 
variable  image).  However,  these  azimuth  ambiguities  were  distinguished  from  the  actual  targets, 
as  shown  in  Figure  20  (the  HH  image),  which  zooms  in  on  the  areas  of  interest.  The  detection 
results  are  the  red  squares,  the  targets  are  green  rectangles  and  the  ambiguities  are  yellow 
rectangles.  These  initial  results  demonstrate  that  small  ships  are  more  likely  to  be  detected  using 
the  KCFAR  Quad  detector. 
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□  KCFARQuad 


AIS 


Figure  19:  KCFAR  detection  results  of  the  Strait  of  Gibraltar  in  a  geodetic  projection. 

Original  RADARS AT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 
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KCFARQuad 

1 — ! - 1  target 

1  _ 1  Azimuth  ambiguity 

AIS 

Figure  20:  KCFAR  detection  results,  a)  and  b)  land  ambiguities;  c)  a  target  and  its  ambiguity. 

Original  RADARSAT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd.,  2008 
-  All  Rights  Reserved. 
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9  Conclusions 


This  study  has  evaluated  the  baseline  performance  of  ship  detection  algorithms  applied  to 
RADARS AT-2  FQ  data.  The  algorithms  were  originally  developed  using  airborne  PolSAR 
images  from  the  CV-580  system  and  this  study  verifies  that  the  same  algorithms,  which  are  based 
on  the  likelihood  ratio  test,  are  applicable  to  RADARS  AT-2  images  as  well. 

Ship  detection  performance  for  various  SAR  systems,  including  quad  polarization,  HH  +  F1V  with 
amplitude  only,  and  single  polarization  HH  has  been  estimated.  In  addition,  a  CP,  right  circular 
transmit  and  linear  receive,  SAR  system  was  simulated  using  measured  data  from  RADARS  AT-2 
FQ  Mode.  The  polarimetric  SAR  ship  detection  algorithms  were  also  applied  to  this  CP  system. 

It  has  been  shown  that  when  a  CP  system  is  simulated  using  RADARSAT-2  FQ  data,  the  noise 
power  in  each  CP  channel  is  twice  the  noise  power  in  each  RADARSAT-2  channel.  However, 
the  two  channels  of  the  CP  system  are  dominated  by  the  co-polarization  channels  of  the 
RADARSAT-2  image.  Since  the  ocean  clutter  of  the  co-polarization  channels  is  about  30  dB 
above  the  noise  floor,  the  effective  target  to  clutter  ratio  will  be  affected  only  slightly  by  the 
thermal  noise  floor.  The  system  performance  is  approximately  determined  by  the  effective  target 
to  clutter  ratio  of  the  signal  components.  Therefore,  the  effect  of  the  RADARSAT-2  system  noise 
on  the  simulated  CP  system  detection  performance  can  be  ignored. 

The  quad  polarization  system  provides  the  best  detection  in  this  application;  the  dual  polarization 
systems,  simulated  CP  and  HH  +  HV  with  amplitude  only,  provide  better  detection  than  single 
channel  systems.  For  five  of  the  test  targets  studied,  the  CP  SAR  system  provides  significantly 
better  detection  than  the  HH  +  HV  amplitude  only  system.  For  the  other  three  cases,  the  HH  + 
VH  system  provides  slightly  better  detection  than  the  CP  system.  It  is  judged  that  the  CP  system 
is  better  than  the  HH  +  HV  amplitude  only  system.  The  CP  SAR  system  therefore  has  potential 
for  improved  ship  detection.  At  the  same  time,  a  dual  polarization  system  will  have  wider  swath 
coverage  than  that  of  a  quad  polarization  system  having  the  same  resolution. 

It  was  also  verified  that  the  phase  difference  between  the  HV  and  the  VH  channels  can  aid  in 
distinguishing  targets  and  azimuth  ambiguities,  consequently,  reducing  the  false  alarm  rate. 

The  new  KCFAR  ship  detector  (KCFAR  Quad)  based  on  feeding  the  PolSAR  decision  variable 
derived  from  the  likelihood  ratio  test  into  a  KCFAR  detector,  has  been  tested.  The  algorithms  for 
calculating  the  decision  variable  from  quad  polarization  data  was  implemented  as  the  required 
parameter  of  the  KCFAR  detector  in  ADSS.  The  initial  results  demonstrate  the  applications  of 
this  detector  for  ship  detection  with  quad  polarization  data. 

The  RADARSAT  Constellation  Mission  will  provide  CP  SAR  as  an  operational  mode.  It  is 
recommended  that  this  mode  be  considered  for  wide  area  surveillance  since  it  will  provide  wider 
swath  coverage  than  quad-pol  mode,  yet  it  has  additional  phase  information  which  conventional 
dual  polarization  does  not  have.  The  phase  information  will  aid  in  improving  ship  detection  to 
reduce  the  probability  of  false  alarms. 
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Annex  A  Detection  performance 


The  detection  performance  of  data  obtained  from  RADARSAT-2  FQ  Mode  as  well  as  the 
simulated  CL  data  have  been  studied,  and  the  ROC  curves  determined.  The  detection 
performance  analysis  has  been  discussed  in  Section  7.5  and  the  comparison  of  the  performance  of 
various  systems  has  been  listed  in  Table  4.  This  Annex  presents  additional  ROC  curves.  Figures 
a)  present  the  ROC  curves,  Figures  b)  are  the  images  of  the  |HH|  channel  and  Figures  c)  are  the 
images  of  the  |X1|  component  of  the  simulated  CL  data.  The  results  shown  in  this  annex  and  in 
Section  7.5  clearly  demonstrate  that  a  quad  polarization  system  (triangles)  provides  the  best 
detection  for  the  cases  studied.  A  dual  polarization  system,  the  F1F1+HV  amplitude  only  system 
(circles)  and  the  simulated  CP  data  (squares)  provide  better  detection  than  a  system  operating  in  a 
single  HH  polarization  (asterisks).  For  five  of  eight  cases  studied,  the  simulated  CL  data  provides 
better  performance  than  the  F1F1+F1V  amplitude  only  system.  For  other  three  cases,  the  F1F1+F1V 
with  amplitude  only  system  provides  slightly  better  performance  than  the  simulated  CL  data. 
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Figure  21:  a)  Detection  performance  for  ship  1;  b)  \HH\  image  and  c)  simulated  \X/\  image. 

Original  RADARSAT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 
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Figure  22:  a)  Detection  performance  for  ship  4;  b)  \HH\  image  and  c)  simulated  \Xj\  image. 

Original  RADARS AT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 
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Figure  23:  a)  Detection  performance  for  ship  7;  b)  \HH\  image  and  c)  simulated  \X/\  image. 

Original  RADARSAT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 
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a) 


Figure  24:  a)  Detection  performance  for  ship  8;  b)  \HH\  image  and  c)  simulated  \X/\  image. 


Original  RADARSAT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 
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Figure  25:  a)  Detection  performance  for  ship  9;  b)  \HH\  image  and  c)  simulated  \X/\  image. 

Original  RADARSAT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 
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Figure  26:  a)  Detection  performance  for  ship  10;  b)  \HH\  image  and  c)  simulated  \Xt\ 
image. 

Original  RADARSAT-2  Data  and  Products  ©MacDonald,  Dettwiler  and  Associates  Ltd., 
2008  -  All  Rights  Reserved. 
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